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Abstract. Spatial birth and death processes are obtained as solutions of a sys- 
tem of stochastic equations. The processes are required to be locally finite, but 
may involve an infinite population over the full (noncompact) type space. Con- 
ditions are given for existence and uniqueness of such solutions, and for temporal 
and spatial ergodicity. For birth and death processes with constant death rate, a 
sub-criticality condition on the birth rate implies that the process is ergodic and 
converges exponentially fast to the stationary distribution. 



1. Introduction 

Spatial birth and death processes in which the birth and death rates depend on 
the configuration of the system were first studied by Preston (1975). His approach 
was to consider the solution of the backward Kolmogorov equation, and he worked 
under the restriction that there were only a finite number of individuals alive at 
any time. Under certain conditions, the processes exist and are temporally ergodic, 
that is, there exists a unique stationary distribution. The more general setting 
considered here requires only that the number of points alive in any compact set 
remains finite at all times. 

Specifically, we assume that our population is represented as a countable subset 
of points in a complete, separable metric space S (typically, S C R d ). We will 
identify the subset with the counting measure rj given by assigning unit mass to 
each point, that is, r](B) is the number of points in a set B G B(S). (B(S) will 
denote the Borel subsets of S.) We will use the terms point process and random 
counting measure interchangeably. With this identification in mind, let Af(S) be 
the collection of counting measures on the metric space S. The state space for 
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our process will be some subset of f\f(S). All processes and random variables are 
defined on a complete probability space (O, J 7 , P). 

The spatial birth and death process is specified in terms of non-negative functions 
X : S x N(S) — > [0, oo ) and 6 : S x Af(S) — > [0, oo) and a reference measure /? on 
5 (typically Lebesgue measure ma, if 5 C A is the birth rate and (5 the death 
rate. If the point configuration at time t is r) G J\f(S), then the probability that 
a point in a set B C 5 is added to the configuration in the next time interval 
of length At is approximately J B X(x, r])p(dx)At and the probability that a point 
x G r\ is deleted from the configuration in the next time interval of length At is 
approximately 5(x,t])At. Under these assumptions, the generator of the process 
should be of the form 

AF( V ) = J (F(r ] +5 X )-F(r ] ))\(x,r 1 )(3(dx) + J (F( n -5 X )-F(ri))5(x,ri)r)(dx) (1.1) 

for F in an appropriate domain. 

Following the work of Preston, spatial birth and death processes quickly found 
application in statistics when Ripley (1977) observed that spatial point patterns 
could be simulated by constructing a spatial birth and death process having the 
distribution of the desired pattern as its stationary distribution and then simulating 
the birth and death process for a long time, a procedure now known as Markov chain 
Monte Carlo. 

The two best-known classes of spatial point processes are Poisson random mea- 
sures and Gibbs distributions. 



1.1. Poisson random measures. Let (3 be a cr-finite measure on S, (S,ds) a com- 
plete, separable metric space. £ is a Poisson random measure on S with mean 
measure (3 if for each B G B(S), £(B) has a Poisson distribution with expectation 
f3(B) and £(B) and £(C) arc independent if B n C = 0. Taking A = 5 = 1, then 
the Poisson random measure with mean measure (3 gives the unique stationary 
distribution for the birth and death process with generator 



= J (F( v + 5 X ) - F( V ))(3(dx) + J 



AF( V ) - / (F( V + S x ) - F( v ))f3(dx) + / (F^ - 5 X ) - F( V ))r)(dx). (1.2) 



Letting ^ denote this distribution, the stationarity can be checked by verifying 
that 



/ 



AF(r,)$(dr,) = 0. 
This assertion follows from the standard identity 

E[f h(Z-5 x ,x)Z(dx)]=E[[ h(t,x)(3(dx)}. (1.3) 
See Daley and Vere- Jones (1988), p. 188, Equation (6.4.11). 

1.2. Gibbs distributions. Assume that (3(S) < oo. Consider the class of spatial point 
processes specified through a density (Radon-Nikodym derivative) with respect to 
a Poisson point process with mean measure (3, that is, the distribution of the point 
process is given by 

^H(dv) = ^—e- H M^(d V ), (1.4) 
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where H(rj) is referred to as the energy function, Zp t u is a normalizing constant, 
and pPp is the law of a Poisson process with mean measure /?. Therefore, the state 
space for this process is S — {n £ Af(S); H(rj) < oo}, the set of configurations with 
positive density. We assume that H is hereditary in the sense of Ripley (1977), 
that is H(rf) < oo and r] C r] implies H(fj) < oo. Ripley showed that such a 
measure Pfj.H is the stationary distribution of a spatial birth and death process. In 
fact, there is more than one birth and death process that has /U/3,i? as a stationary 
distribution; we simply require that \{x,n) > if H(rj + 8 X ) < oo and that A and 
5 satisfy 

\{x, n)e~ H ^ = S(x, n + 5 x )e- H(r >+ s *K (1.5) 

This equation is a detailed balance condition which ensures that births from n to 
1] + S x match deaths from rj + 5 X to r\ and that the process is time-reversible with 
(1.4) as its stationary distribution. Again, this assertion can be verified by showing 
that 

This identity again follows from (1.3). 

Notice that equation (1.5) says that any pair of birth and death rates such that 

j^^.^-H^ «.) + «(„)> 

will give rise to a process with stationary distribution given by (1.4). We can always 
take 5(x,rj) — 1, that is, whenever a point is added to the configuration, it lives an 
exponential length of time independently of the configuration of the process. 

For example, consider a spatial point process on a compact set Sctf given by 
a Gibbs distribution with pairwise interaction potential p{x\,X2) > 0, that is, for 

Em c 

H p(v) = ^P{xi,Xj) (1.6) 

i<j 

= \[J J p{x,y)rj(dx)ri(dy) - J p(x,x)rj(dx)}} 

and the distribution of the point process is absolutely continuous with respect to 
the spatially homogeneous Poisson process with constant intensity 1 (or cquiva- 
lently with Lebesgue mean measure) on S. Taking 8(x,rj) = 1 and X(x,rj) — 
exp{— J p(x, y)n(dy)}, the distribution determined by (1.6) is the stationary distri- 
bution for the birth and death process with infinitesimal generator 

AF(rj) = J e- l ^y^ d y\F(r 1 + S x )-F(7 1 ))dx + J (F{ V -8 x )-F(r,))r,(dx). (1.7) 

Another example is the area-interaction point process introduced by Baddeley 
and Van Lieshout (1995). This point process is absolutely continuous with respect 
to the spatial Poisson process with Lebesgue mean measure on S C M. d and 
H(i]) = i](S) log/? — m d (i] (B G), so the Radon-Nikodym derivative is given by 

L(v) - 4p" (S) 7" md(, ' eG) - (1-8) 

Again, Z is the normalizing constant, p and 7 are positive parameters, and G is a 
compact (typically convex) subset of R d referred to as the grain. The set 77 © G is 
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given by 

77 © G = U{x © G;x e r]}. 

The parameter 7 controls the area-interaction among the points of r\: the process 
is attractive if 7 > 1 and repulsive otherwise. (See Lemma 3.6.) If 7 = 1 the 
point process is just the Poisson random measure with mean measure prrid- The 
case 7 > 1 is related to the Widow-Rowlinson model introduced by Widow and 
Rowlinson (1970). The case of area- exclusion corresponds to a suitable limit 7^0. 
A birth and death process with stationary distribution given by the area-interaction 
distribution can be obtained by taking the unit death rate and the birth rate given 

by 

1.3. Overview. The spatial birth and death processes that correspond to the Gibbs 
distributions discussed above involve finite configurations and it is straightforward 
to see that they are uniquely characterized by their birth and death rates, for 
example, as solutions of the martingale problem associated with the generator A 
given in (1.1); however, if the configurations are infinite and the total birth and 
death rates are infinite, the existence and uniqueness of the processes are not so 
clear. In Section 2, we represent these processes as solutions of a system of stochastic 
equations and give conditions for existence and uniqueness of solutions for the 
equations as well as for the corresponding martingale problems. These equations 
are very useful in studying the asymptotic behavior of the birth and death processes, 
including temporal and/or spatial-ergodicity and the speed of convergence to the 
stationary distribution. 

The uniqueness conditions given here are direct analogs of Liggett's (1972) con- 
ditions for existence and uniqueness for lattice indexed interacting particle systems. 
Stochastic equations for lattice indexed systems were formulated in Kurtz (1980) 
using time-changed Poisson processes and existence and uniqueness given under 
Liggett's conditions. Stochastic equations for spatial birth and death processes of 
the type considered here were formulated in Garcia (1995) using a spatial version 
of the time-change approach. Existence and uniqueness were again given under 
analogs of Liggett's conditions. 

One disadvantage to the time-change approach taken in Kurtz (1980) and Garcia 
(1995) is that the filtration to which the process is adapted depends on the solution. 
A stochastic equation that avoids this difficulty can be formulated by representing 
the birth process as a thinning of a Poisson random measure. Intuitively, the 
approach is analogous to the rejection method for simulating random variables. 
The fact that counting processes and more general marked counting processes can 
be obtained by thinning Poisson random measures is well known, particularly in the 
context of simulation. (See, for example, Daley and Vere-Jones (2003), Section 7.5.) 
Stochastic equations exploiting this approach were formulated for lattice systems 
in Kurtz and Protter (1996) and for general spatial birth processes by Massoulie 
(1998). In both cases, uniqueness was obtained under conditions analogous to 
Liggett's. 

Section 3 considers temporal ergodicity for birth and death processes in non- 
compact S (more precisely, S and (3 with (3(S) = 00) and spatial ergodicity for 
S = R d and translation invariant birth rates. We give conditions for ergodicity and 
exponential convergence to the stationary distribution. It is well known that these 
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processes are temporally and spatially ergodic if the birth and death rates are con- 
stant (the stationary measure being Poisson). More generally, in Theorem 3.10, we 
show that if the birth rate satisfies the conditions of Theorem 2.13 with M < 1 and 
the death rate is constant (S = 1), the system is temporally ergodic and for every 
initial distribution, the distribution of the solution converges at an exponential rate 
to the stationary distribution. For S = R d and A translation invariant, from the 
stochastic equation, we see that spatial ergodicity of the initial distribution implies 
spatial ergodicity of the solution at each time < t < oo. Unfortunately, it is not 
clear, in general, how to carry this conclusion over to t — oo, that is, to the limiting 
distribution of the solution, although in the case M < 1, spatial ergodicity holds 
for the unique stationary distribution as well. We give some additional conditions 
under which spatial ergodicity of the limiting distribution can be obtained. 

Fernandez, Ferrari and Garcia (2002) study ergodicity of spatial birth and death 
processes using a graphical representation to construct the stationary distribution 
that is closely related to the stochastic equations we consider here. They give 
conditions for an exponential rate of convergence to the stationary distribution and 
for spatial ergodicity of the stationary distribution similar to those given here, but 
for a more restricted class of models. 

Throughout, C(S) will denote the space of bounded continuous functions on S 
and B(S) the Borel subsets of S. 

The stochastic equations we consider will be driven by a Poisson random measure 
N on U x [0, oo) for an appropriate space U, having mean measure of the form vxnii, 
where mi is Lebesgue measure on [0, oo). Then for B 6 B(U) with f(B) < oo, 
N(B,t) = N(B x [0,t]) is just an ordinary Poisson process with intensity v{B). 
For a filtration {Ft}> we sa Y that N is compatible with {J-t} if and only if for each 
B e B(U) with v(B) < oo, N(B, •) is {:F t }-adapted and N(B,t + s) - N(B, t) is 
independent of T t for s, t > 0. 

2. Spatial birth and death processes as solutions of stochastic equations 

A birth and death process as described in the previous section can be represented 
as the solution of a system of stochastic equations. The approach is similar to Garcia 
(1995) where such processes were obtained as solutions of time-change equations. 
We assume that the individuals in the birth and death process are represented by 
points in a Polish space S. Typically, S will be M. d , Z d , or a subset of one of these, 
but we do not rule out more general spaces. Let K\ C K<i C • • • satisfy UkKk = S, 
and let Ck G C(S) satisfy Ck > and mf xe K k Ck{%) > 0. Af(S) will denote the 
collection of counting measures on S and S will denote {£ £ Af{S) : J s Ck(x)C t (dx) < 
oo,fc=l,2,...}. Without loss of generality, we can assume that c\ < C2 < • • • . Let 
C = {/ e C{S) : |/| < acfe for some k and a > 0}, and topologizc S by the weak* 
topology generated by C, that is, Cn —* C if an d only if J s fd( n — > J s fd( for all 
/ G C. (Note that C is linear and that, with this topology, S is Polish.) D$[0, oo) 
will denote the space of cadlag S- valued functions with the Skorohod (Ji) topology. 
We assume that A and 5 are nonnegative, Borel measurable functions on S x C. 

Let P be a cr-finite, Borel measure on S. We assume 

Condition 2.1. For each compact /C C S, the birth rate A satisfies 




t>0, k = l,2, ... 



(2.1) 
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and 

5{x,() <oo, (eS, xe(. (2.2) 
We also assume that A and S satisfy the following continuity condition. 
Condition 2.2. If 

lim f c k (x)\t n -t\(dx)=0, (2.3) 

for each k = 1, 2, . . ., then 

X(x,() = lim X(x,C n ), 6(x,0= Urn 5(x,C n ). (2.4) 

n^oo n^oo 

Note that since (2.3) implies £„ converges to C in S, the continuity condition 
(2.4) is weaker than continuity in S: however, we have the following condition 
under which convergence in S implies (2.3). 

Lemma 2.3. Suppose Co, (1,(2, ■ ■ ■ ^ <S and Cn < Co? n = 1, 2, . . .. If ( n — > C m <S, 
tften Zio/ds. 

Proof. ( n < Co implies that, considered as a measure, Cn « Co an d ^ = 
^TrfjT < 1, almost everywhere Co- Furthermore, Cn — > C m S implies ( n ({x}) — > 
C({a;}) for each x G Co, since the support of Co consists of a countable collection of 
isolated points. Consequently, 

C ^)>c fc (x)|^-if»H0, 

Co(M) Co(M) 

and since J s c k (x)(o(dx) < 00, the dominated convergence theorem implies 

lim / c fc (x)|Cn-CI(^)= I™ / c k (x)\^^ - ^j^r\C (dx) = 0. 
n ^°°Js n ^ooJ s Co(W) Ca({x}) 



□ 



Lemma 2.4. Suppose H C 5 and Co S 5 satisfy ( < (0, ( E Tt. If S satisfies (2.2) 
and Condition 2.2, then 

sup (5(x, C) < 00. 
(en 

Proof. Tt is relatively compact in S, so any sequence {Cn} C H has a subsequence 
that converges in S and, by Lemma 2.3, satisfies (2.3). Fix x G S, and let {Cn} sat- 
isfy lim„^oo S(x, Cn) = sup^ eW S(x, C). Then there is a subsequence that converges 
to some C G S and hence, sup^ gH S(x, () = S(x, () < 00. □ 

Lemma 2.5. Suppose that for each i£S, i/iere exists k(x) such that X(x, ( + 6 y ) = 
X(x, C) for y £ K k ( x y Then X satisfies Condition 2.2 and similarly for 5. 

Proof. Note that (65 implies C{Kk{x)) < 00 and that Cn — * C implies that for 
n sufficiently large, Cn restricted to K k ( x ) coincides with C restricted to K k ( x ), so 
A(x,C„) = A(x,C). □ 

Let TV be a Poisson random measure on S x [0, oo) 3 with mean measure (3{dx) x 
ds x e~ r dr x du. Let rjo be an <S-valued random variable independent of N, and let 
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fj be the point process on S x [0, oo) obtained by associating to each "count" in 770 
an independent, unit exponential random variable, that is, for 770 = $xi> set 

oo 

rjo =Y, 6 (xi,n), (2-5) 

i=l 

where the {n} are independent unit exponentials, independent of 770 and N. The 
birth and death process 77 should satisfy a stochastic equation of the form 

Vt(B) = / %>,A(a:,?7 S _)] r * 5(x,ri v ) dv,oo) {r)N(dx, ds, dr, du) 

JSx[0,t]x[0,oo) 2 Js 

+ / \p n 5(x,r, s )ds,oo){ r )rio{dx,dr). (2.6) 

JBx[0,oo) 

To be precise, let 77 be a process with sample paths in Dg[0, 00) that is adapted to 
a filtration {T t } with respect to which N is compatible. (Note that (2.1) ensures 
that the integral with respect to TV on the right exists and determines an 5-valued 
random variable, and the continuity condition (2.4) and the finiteness of 5(x,Q 
ensure that S(s,r] t ) is a cadlag function of t, so that the integrals J S(x,r] v )dv 
exist.) Then 77 is a solution of (2.6) if and only if the identity (2.6) holds almost 
surely for all B 6 B(S) and t > (allowing 00 = 00). 

Lemma 2.6. Suppose Condition 2.1 holds. If t] is a solution of (2.6), then for 
each T > 0, 

/ / Ck(x)\(x, 7] s )P(dx)ds < 00 a.s., (2.7) 
Jo Js 

77^ defined by 

Vt( B ) = / l{o.\( x .- ns -)](u)N(dx,ds,dr,du) 

JBx[0,T]x[0,oo) 2 

is an element of S, 

Vt<Vr + Va, 0<t<T, 

and 



lim / cfc(x) \r] s - rjt \ (dx) = 0, t > 0. 
s ^t+ Js 



Proof. Since for almost every ui e Q, the closure of {77^ : < s < T} is compact, 
Condition 2.1 implies 



is 

and hence (2.7). Letting 



sup / Ck(x)X(x, n s )(3(dx) < 00 a. 
8<tJs 



in{{t : / / c k (x)\(x, n s )f3(dx)ds > c}, 
Jo Js 



we have 

rT/\T c r 

E[ c k (x)r]^ ATc (dx)} = E[ I I c k (x)X(x, r ls )(3(dx)ds] < c, 
Js Jo Js 

and since lim^oc t c = 00 a.s., it follows that f s Ck{x)i]^(dx) < 00, a.s. implying 
77J e S a.s. The last statement then follows by Lemma 2.3. □ 
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If r\ is a solution of (2.6) and a point at x was born at time s < t, then the 
"residual clock time" r — f* 5(x,r] v )dv is an ^-measurable random variable. In 
particular, the counting-measure-valued process given by 

rj t (B x D) = I \o,\(x, v .-)\{u)^-D{r - I 5(x,7] v ) dv)N(dx 1 ds,dr 1 du) 

JBx[l),(]x[0,oc) 2 Js 

+ I l D {r- ( 5{x,ri s -)ds)rjQ{dx,dr) (2.8) 

JBx[0,oo) JO 

is {^tj-adapted. 

Let S denote the collection of counting measures ( on S x [0, oo) such that 
((• x [0, oo)) G S. We can formulate an alternative equation for the 5-valued 
process r} by requiring that 



f(x,r)rjt(dx,dr) (2.9) 

Sx[0,oc) 

/ f{x,r)r} {dx,dr) 

JSx[0,oc) 

f(x, r)l[o,\( x . Vs _)] (u)N(dx, ds, dr, du) 

l,oo) 2 

/'/ 

JO JSx[0,< 



Sx[0,t]x[0,oo) 2 

5(x, T] s )f r (x, r)rj s (dx, dr)ds, 



!,Oo) 



for all / G C, where C is the collection of / G C(S x [0, oo)) such that f r = J^/ G 
C(S x [0, oo)), f(x,0) = 0, sup r /(•, r)|, sup r \ f r (-,r)\ G C, and there exists r/ > 
such that / r (x, r) = for r > r/. Note that if / G C and 



"(x,r) = / |/ r (a;,u)|dt 
Jo 



(2.10) 



then /* G C. In (2.9), % can be any S- valued random variable that is independent 
of N. 



2.1. Martingale problems. Let V{A) be the collection of functions F of the form 
F(() = e _ ^ Sx [°.°°) f( x > r )(( dx > dr ) j f or non-negative / G C. Suppose that ?/ is a solution 
of (2.9) with sample paths in f^[0, oo). Assuming Condition 2.1, 

/ / c k (x)X(x,r, s )/3(dx)ds < oo, fc = 1,2,.... (2.11) 
Jo Js 
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By Ito's formula 



F(%) = F(r} ) + f F(ff a -)(e- f ^ - l)i { o Mx , v ..)]{u)N{dx,d8,dr,du) 

J Sx[0,t]x[0,oo) 2 



+ / F(Vs) f 5{x, Vs )f r (x,r)rj s (dx,dr)ds, (2.12) 

JO Jsx[0,oo) 



= F(rjo) + f F(rf s -)(e- f ^ r) - l^^.^Nidx^d^du) 

JSx[0,t]x[0,oo) 2 

+ I F{rj s )( [ X(x, Vs )(e-f (x ^ -l)e- r l3(dx)dr 

JO V JSx[0,oo) 

+ / S(x,r] s )f r (x,r)r} s (dx,dr))ds, 

JSx[0,co) ' 

where N(dx,ds,dr,du) = N(dx,ds,dr,du) — (3(dx) x ds x e~ r dr x du. It follows 
from (2.11) that the stochastic integral term on the right is a local martingale, and 
since /* given by (2.10) is in C, it follows that 



/7 

Jo JSx[0. 



S(x, rj s )\f r (x, r)\ff s (dx, dr)ds < oo, t > 0. (2.13) 

oo) 



Consequently, defining 



AF(C) = F(C)( \{xX){e- f[x > r) - l)e' r I3{dx)dr 

W,Sx[0,oo) 



f 6(x,Of r (x,r)((dx,dr)), (2.14) 

JsxlO.oc) ' 



'Sx[0,oo) 

any solution of (2.9) must be a solution of the local martingale problem for A. 
We say that rj is a solution of the local martingale problem for A if there exists a 
filtration {Ft} such that rj is {^j-adapted and 

M F (t) = F(fj t ) - F(rj ) - f AF(rj s )ds (2.15) 

Jo 

is a {^-"tj-local martingale for each F e D(A), that is, for each F of the form 
F(() = e~ / f EC, f > 0. In particular, let J(x) = sup r f(x, r) and 

t/ jC = inf{t : / / f(x)\(x,r) s )l3(dx)ds > c}. 
Jo Js 

Then Mp(- A r/ jC ) is a martingale. Note that r/ jC is a {^j-stopping time. 

Conversely, if fj is a solution of the local martingale problem for A with sam- 
ple paths in Dg[0,oo), then under Condition 2.1, (2.11) and (2.13) hold. If 7 6 
Cr[0, 00) has compact support and f(x,r) = j(u)duck(x), then / € C and it 
follows that 

/ / Cfe (x)S(x, rj s )fj s (dx, dr)ds < 00, i > 0. 

To formulate the main theorem of this section, we need to introduce the notion of 
a weak solution of a stochastic equation. 
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Definition 2.7. A stochastic process f\ with sample paths in D§[0, oo) is a weak 
solution of (2.8) if there exists a probability space (O,^ 7 , P), a Poisson random 
measure ./V on Sx [0, oo) 3 with mean measure [3(dx)xdsxe~ r drxdu and a stochastic 
process 77 defined on (Q,T,P), such that 77 and 77 have the same distribution on 
Dg[0, 00), rj is adapted to a filtration with respect to which N is compatible, and 
N and 77 satisfy (2.8). 

Theorem 2.8. Suppose that A and 5 satisfy Conditions 2.1 and 2.2. Then each 
solution of the stochastic equation (2.8) (or equivalently, (2.9)) is a solution of the 
local martingale problem for A defined by (2.14), and each solution of the local 
martingale problem for A is a weak solution of the stochastic equation. 

Proof. The first part of the theorem follows from the discussion above. 

To prove the second part, we apply a Markov mapping result of Kurtz (1998). 
Let {Di} C B(S x [0, oo) 2 ) be countable, closed under intersections, generate B(S x 
[0, oo) 2 ), and satisfy J D (3{dx)e~ r drds < 00. Then N is completely determined by 
N(Di,t). Define 

Zi(t) = Z i (0)(-l) JV < 2, " t \ 
where Z»(0) is ±1. Note that 

1 r* 

N(D h t) = -- Zi(a-)dZi(8), (2.16) 

and if the Z,(0) arc iid with P{Z t (0) = 1} = P{Z t (Q) = -1} = \ and independent 
of N, then for each t > 0, the Zi(t) are iid and independent of iV. For z G { — 1, 1} 00 , 
we will let denote 

^_ 1 ^ D (x,r,u) z = ((_i)i Dl (^,«) 0l5 (-l) 1 ^^-'-.")^, . . .). 
Then Z = [Z\, Z 2 , . . .) is a solution of the martingale problems for 

CG(z)= I (G((-l) lD(x < r > u) z)-G(z))e- r drdul3(dx). 

We can take the domain for C to be the collection of functions that depend on only 
finitely many coordinates of z. With this domain, the martingale problem for C is 
well-posed. 

If 7/ is a solution of (2.9), then (77, Z) is a solution of the local martingale problem 

for 

k{FG){lz) 

= F(()( [ ((l [0 , A(x , C)] (77) e -^^'-)+l (A(;c , C) , oo) (7 i ))G((-l) 1 -^)z) 

-G(z)y- r (3(dx)drdu 
-G(z) [ 5(x,0fr(x,r){(dx,drj). (2.17) 

JSx[0,oo) ' 
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Let fj be a solution of the local martingale problem for A. For a = (a\, az, . . .) 
with a k > 0, k = 1, 2, . . ., define 

7V,(t) = 



°° r f- 

inf {m : / lvVflt / c k (x)X(x, fj 8 )(3(dx) 4 
fe=1 Us 

OO r „ 

H„(C) = 1 V V a k / c k (x)X(x, C)p(dx) 
Us 



c k {x)S(x,fj s )fj s (dx) 



SxS 



ds > t}, 



c k {x)8(x,()((dx) 



SxS 



(2.18) 



and 77" = ?? TQ (t). Then fj a is a solution of the martingale problem for 

For F <E ©(^4), -4 a -F is bounded, and wc can select a n = (a™,^, . . .) so that 
a" > a™ +1 and T a » (t) -» f a.s. 

Let /j(rfz) = nfeLi(|^{-i}(^ z fc) + ^{l}^*))) an d set °G = J Gdfi. Then 



/ A(FG)(C,z)m(^) = c g AF(C), 



and more generally, 



AA(FG)(C,z)n(dz) 



c G —AF{Q. 
n a 



Applying Corollary 3.5 of Kurtz (1998) to H~ 1 A for each a, we conclude that if 
fj is a solution of the local martingale problem for A, then there exists a solution 
(fj, Z) of the local martingale problem for A such that fj and fj have the same 
distribution. Finally, applying (2.16), wc can construct the corresponding Poisson 
random measure TV and show that fj and N satisfy (2.9). □ 

The natural (local) martingale problem for r\ is really the martingale problem for 
A given by (1.1); however, there will be solutions fj of the local martingale problem 
for A (and hence of the stochastic equation) such that the corresponding 77 is not 
a solution of the local martingale problem for A. Intuitively, conditioned on = 
a{j] s : s < t), the residual clock times should be independent unit exponentials, 
independent of J^. That need not be the case, since we are free to pick the residual 
clock times at time zero in any way we please. It also need not be the case if the 
solution of the martingale problem fails to be unique. The following results clarify 
the relationship between the martingale problems for A and A. 

Proposition 2.9. Suppose that A and S satisfy Conditions 2.1 and 2.2. If fj is 
a solution of the local martingale problem for A and at each time t, the residual 
clock times are independent of J-^ and are independent unit exponentials, then rj is 
a solution of the local martingale problem for A. 

Proof. By assumption, wc can write fj t = ^2i^(Xi(t),Ri(t))i where the Ri(t) are 
independent unit exponentials, independent of ', and in particular, independent 



of n t . For / e C and F(Q 



-J. 



Sx[0,oo) 



f(x,r)C(dx,dr) 



since (2.15) can be localized 



by {.F^j-stopping times, it follows that 

E[F{%)\^] - E[F(ffo\^} f E[AF( Vs )\^}ds 

Jo 
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is a {jF^j-local martingale. By the independence of the Ri(t), 

r-OO 

E[F{r} t )\^} = II / e~ nXi ^^e- r dr = e~ *s fl(*)«J«(<fa) = 
i J o 

where g is defined so that e^ 9 ^ = f °° e~^ x ' r ^e~ r dr. Integrating by parts gives 

/>OC />oo 

/ e- f(x > r) f r (x,r)e- r dr = 1 - / e- f{ - x ^e~ r dr = 1 - e~^ x \ 
Jo Jo 



and hence 



E[AF( Vs )\rj] = G(rj s ) f X(x, Vs )(e-^-l)e- r p(dx)dr 

JSxlO.oo) 



'Sx[0,oo) 

/>oc 

e -/W(-),r) e -r dr 




/>OC 

/ e -^ x ^S(X j (s),r ls )f r (X j (s),r)e- r dr 
Jo 

G{ Va ) I X(x, Vs )(e- g{x) -l)0(dx) 




■ e - 9 W(-)) I 5(X,-(«),f, a )(l-e-»WW)) 

and the proposition follows. □ 
We have the following converse for the previous proposition. 

Theorem 2.10. Suppose that A and S satisfy Conditions 2.1 and 2.2. If r] is a 
solution of the local martingale problem for A, then there exists a solution rj of the 
local martingale problem for A such that rj and rj(- x [0, oo)) have the same distribu- 
tion on Ds[0, oo) and at each time t > 0, the residual clock times are independent, 
unit exponentials that are independent of ■ 

Proof. For ( = S Xi G S, let a((, d() denote the distribution on S of J^i fyxi,T<)> 
where the n are independent, unit exponential random variables. Then, by the 
calculation in the proof of Proposition 2.9, 

G(C) = / F(0a(Cd() AG(C) = f AF(C)a(C,dC), 

for F G V{A). More generally, A a G(() = J § A a F(()a{(,d(), where A a is defined 
as in (2.18). The theorem then follows by Corollary 3.5 of Kurtz (1998). □ 

Corollary 2.11. Let v G V{S), and define v G V{S) by 



hdv= \ / /i(C)a(C,dCMO- 

5 Js J§ 

If uniqueness holds for the martingale problem for {A,v), or equivalently, weak 
uniqueness holds for the stochastic equation (2.9), then uniqueness holds for the 
martingale problem for (A, v) . 



12 



Nancy L. Garcia and Thomas G. Kurtz 



Proof. If rj is a solution of the martingale problem for (A, v), then Theorem 2.10 
gives a corresponding solution of the martingale problem for (A,v). Uniqueness for 
the latter then implies uniqueness of the former. □ 



2.2. Existence. We now turn to the question of existence of solutions of (2.6). We 
assume that Conditions 2.1 and 2.2 hold. The pair (A, S) will be called attractive if 
Ci C C2 implies \(x,(i) < C2) and 5(x,(i) > 5(x, £2)- If (A, 8) is attractive and 
we set rp = 0, then if 1 defined by 

% +1 (B) = / lo.ACc.ij?-)] (")%**(* , v s)dv t00) (r)N(dx,d8,dr,du) 

iBx[l),i]x[0,x) 2 

«(*,„») ds,oc)( r )vo(dx, dr) (2.19) 



'Bx[0,oo) 

is monotone increasing and either 77™ converges to a process with values in S, or 

/ [ c k {x)\{x,r%)(3{dx)ds 00, (2.20) 
Jo Js 

for some T and fc. To see this, let 

t c " = inf{t : / / c fe (a;)A(a;,r;™)/3(dx)ds > c}. 
Jo is 

Then 

E[ sup / c fe (a;)^ +1 (dx) - / \f* s(x. n ™) ds.oo)( r )Vo(dx, dr) J] 

KTAt" \^Js JAx[0,oo) Wo ^ ' s ' / 

<E[ / c fc (a;)A(a;,r£)/3(da;)ds] 

Jo Js 

<c, 

and r 1 > r c 2 > • • • . Either 

lim lim r c " = 00 (2.21) 



or (2.20) holds for some T. 

If (2.21) holds almost surely, the limit 77°° is the minimal solution of (2.6) in the 
sense that any other solution 77 will satisfy rj%°(B) < rj t {B) for all B E B(S) and 
t>0. 

For an arbitrary pair (A, 8) satisfying Conditions 2.1 and 2.2, we define an at- 
tractive pair by setting 

X(x, C) = sup X(x, (') 8(x, C) - inf S(x, (')■ 
C'CC 

Let 770 be an 5-valued random variable independent of N, and let fjo be defined as 
in (2.5). We assume that A satisfies (2.1), which implies 



c k (x)\(x,()/3(dx) < 00, (e5,fc = l,2,..., (2.22) 
and that there exists a solution rj for the pair (A, 8). 
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We consider a different sequence of approximate equations. Let {K n } be the 
sets in the definition of C, and let rf 1 satisfy 

Vt(B) = 

/ MH=o,%-nK n n^_)]( u )\j' 5(x,rj v nK n n v z)dv,oo)( r ) N ( dx ' da, dr, du) 

■/Bx[0,t]x[0,oo) 2 Us ' 

+ / 1 {jU(x,r jv nK n nr,n )dVt0o) {r)f}o(dx,dr). (2.23) 

JBx[0,oo) 

Existence and uniqueness for this equation follow from the fact that only finitely 
many births can occur in a bounded time interval in K n . Consequently, the equation 
can be solved from one such birth to the next. Since X(x, r] s _r\K n r\r]™_) < X(x,fJ s _) 
and S(x,rj s n K n n rf£) > 5(x,rj s ), it follows that 77™ C rj t and hence that 

V?(B) = / MaM^,K n nr l ^](u)l(ft S ( x .K n nrr)dv.oo)( r ) N ( dx ^ d ^ dr ^ du ) 

JBx[0,t]x[0,oo) 2 Js 

+ / \fls(x,K n nn2)dv,oo)( r )Vo( dx > dr )- ( 2 - 24 ) 

Also, note that for g e C, 
Us 

g(x)5{x,K n nr 1 l l ) V :{dx)d S < J g{x)r\ Q ^ {x - a )] {u)N(dx,ds,dr,du) < 00. 

(2.25) 

Define F(£) = e" J sx[»,») f{x,r)C(dx,dr) ^ j g £ nonnegative. Setting 

A„F(C) - F(C) ( / A(x, Jf„ n C)(e- /(a: ' r) - lK^daOdr 

V JSx[0,oo) 



as in (2.12), 



+ / S(x,K n nOfr(x,r)((dx,dr)), (2.26) 



/ 8{x, 

JSx[0,oo) 

F{rj?) - F(rf ) - C A n F(r%)ds 
Jo 



'0 

is a local martingale. 

Uniqueness for (2.24) implies that the residual clock times at time t are condi- 
tionally independent, unit exponentials given J 7 ?. Consequently, as in Proposition 
2.9, for G(C) = e-isff^X^), g eC nonnegative, and 

A n G(C) = J (G((+6 x )-G(())\(x,K n nO(3(dx)+ J (G(C-8 x )-G(Q)5(x,K n nQC(dx), 



GOT?) - Gfag) - / A n G{rt n s )ds (2.27) 
Jo 

is a local martingale. Exploiting the fact that 77™ Crj t , the relative compactness of 
{i] n } 7 in the sense of convergence in distribution in Dq[0, 00) follows. 

Proposition 2.12. Suppose that Conditions 2.1 and 2.2 hold. If(x,Q — ► A (a:, C) 
anrf (x, C) — > S(x, Q are continuous on S x C, then Q — > AG(£) is continuous, and 
any limit point of {rj n } is a solution of the local martingale problem for A, and 
hence a weak solution of (2.8). 



14 



Nancy L. Garcia and Thomas G. Kurtz 



Proof. By (2.22), we can select so that 

T(t) = [ lvVat [ c k (x)X(x,rj s )/3(dx)ds < oo, Vt > a.s. 
Jo I Js 

and by (2.25), it follows that 

T m = inf{i : T(t) > m} 

is a localizing sequence for (2.27) for all g and n. The estimates also give the 
necessary uniform integrability to ensure that limit points of (2.27) are local mar- 
tingales. □ 

2.3. Existence and Uniqueness. If sup^ g5 J s X(x, Q[3{dx) < oo, then a solution of 
(2.6) has only finitely many births per unit time and it is easy to see that (2.6) has 
a unique solution. Condition 2.1, however, only ensures that there are finitely many 
births per unit time in each Kk, and uniqueness requires additional conditions. The 
conditions we use are essentially the same as those used for existence and uniqueness 
of the solution of the time change system in Garcia (1995). From now on, we are 
going to assume that 5(x, n) = 1, for all x S S and n 6 S. 

Let N be a Poisson random measure on S x [0, oo) 3 with mean measure (3(dx) x 
ds x e~ r dr x du. Let 770 be an 5-valued random variable independent of N, and let 
rjo be defined as in (2.5). Suppose {T t } is a filtration such that rjo is .^-measurable 
and is {^j-compatible. We consider the equation 



Vt(B) = / li 0i x(x,n 3 -)](u)l{t-s,oo)(r)N(dx,ds,dr,du) 

JBx[0,(]x[0,oc) 2 



/ \ ttao) (r)no(dx,dr). (2.28) 

JBxfO.ool 



'Bx[0,oo) 

Theorem 2.13. Assume Conditions 2.1 and 2.2. Suppose that 

a(x,y) > snp\X(x,n + S y ) - \(x,r))\ 
v 



and that there exists a positive function c such that 

= sup / 

x JS 

Then, there exists a unique solution of (2.28) 



. c(x)a(x, y) n , , . 
M = sup / v ' \ ' Pidy) < 00. 
S c(y) 



Example 2.14. Let d{x,rj) = mi{ds(x, y) : y G n}, where ds is a distance in S 
such that (S,ds) is complete separable metric space. Suppose X(x,rj) = h(d(x,rf)). 
Then a(x,y) = sup r><js / Xi j / \ \h(r) — h(ds(x,y))\. If h is increasing, then a(x,y) = 
h(oo) — h(ds(x,y)) and 

lA^ry 1 ) - X(x, V 2 )\ < J (h(oo) - h(d s (x, y)))^ 1 - V 2 \(dy). 

If h is decreasing, then a(x,y) — h(ds(x,y)) — h(oo) and 

lA^ry 1 ) - X(x, V 2 )\ < J (h(d s (x,y)) - Moo))!*? 1 - V 2 \(dy)- 

Theorem 2.13 is a consequence of the following lemmas that hold under the 
conditions of the theorem. 
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Lemma 2.15. For any 77 1 , ?y 2 £ S we have 

lA^r^-A^ry 2 )! < f a(x,y) - V 2 \(dy). (2.29) 
Js 

Proof. Since 77 1 and ij 2 contain countably many points, there exist {j/i,2/2j ■ • ■} 
and {zi, Z2, ■ ■ ■} such that 

1 J 

i=l j=l 

(where / and J may be infinity) and hence 

\v 1 -V 2 \(B)=J2S m (B)+J2S Zj (B). 
»=i j=i 

By the definition of a and Condition 2.2 

/An jAn 

lACar,^ 1 ) - A(ar,r7 2 )| = lim \X(x, n 1 ) - X(x, n 1 + # y - 5 Z )| (2.30) 

i=l j=l 

(/An JAn 
a(a;, yi) + a(x, Zj) 
»=i j=i 



< 

/s 



/ l^ 1 -rf\{dy). 

Js 



□ 



Define 



Vo(B,t)= \ tjCX>) (r)r}Q(dx,dr). 

JBx[0,oo) 

Let 77 be {^" t }-adapted with sample paths in Ds[0, 00). Then by Condition 2.1 
$Vt(B) = r)o(B,t)+ / l [aMx . Vs)] (u)l (t _ St00) (r)N(dx,ds,dr,du) (2.31) 

JBx[0,t]x[0,oo) 2 

defines a process adapted to {J- t } with sample paths in Dg[0, 00). 

Lemma 2.16. Leti] 1 andrj 2 be adapted to {Ft} and have sample paths in As[0, 00). 
Then 

supc(x)E[ f a(x, y)\W(t) - $?? 2 (*)!(<%)] (2.32) 

x Js 

<M [ sup c(x)E[[ a{x,y)\ril-r£\{dv)]e-<t-')ds. 

Jo x J s 
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Proof. Let C = §rf. Then 

sup c(z)E[[ a{z,x)\^-e t \{dx)\ 

z JS 

< sup c(z)E[ / a(z,a;)|l [0jA(a ,^i )] (u) - l{o,A(s,»7?)](«)l 

z JSx[0,t]x[0,oo) 2 

\t-s,<x)(r)N(dx, ds, dr, du)} 

< supc(z)E[ I a{z,x)\X{x,rr\) - \(x, rr])\e- {t - s) f3{dx)ds] 

z JSx[0,t] 

<supc(z) I a(z,x)E[[ a(x,y)\r]l - 7 1 2 \(dy)}e~ < - t ' s ^(dx)ds 

z JSx[0,t] JS 

< sup c(z) f ^^f3(dx) f sup c(x)E[ [ a(x,y)\r,l-r, 2 s \(dy)}e-^- s Us 

Js c(x) Jo x Js 

< M [ sup c(x)E[ f a(x,y)\ril - ^l(^)]e" ( *" a) <&■ (2.33) 

□ 

Proof. (Theorem 2.13) Uniqueness follows by (2.32) and Gronwall's inequality. 
To prove existence, we proceed by iteration. Let 77° = r] (-, t), and for n > 1, define 

v n+l = Thcn 

sup c(x)E[[ a(x,yM +1 -r,?\(dy)] 

x JS 

<M f S upc(x)E[[ a(x,yM -^Kdyfte-^ds 

JO x Js 

<M 2 f Psupc(x)E[[ aix^M; 1 -^; 2 \(dy)]e-^-^d S2 e-^-^d Sl 

Jo Jo x Js 

<M n f f 1 ... ^ sup c(x)E[f aix^^l^-vUidy) 

Jo Jo Jo x Js 

e -(*»-i-«»)d8 n ... e-^'^ d Sl . 
Therefore, there exists C > such that 

supc(x)E[ a(x,yM +1 ^^\(dy)} < supsupc(x)£;[ / a(x, y)\rr] - V° s \(dy)}, 

x Js n - s<t x Js 

and the convergence of 77™ to a solution of (2.28) follows. □ 



3. Ergodicity for spatial birth and death processes 

3.1. Temporal ergodicity. The statement that a Markov process is ergodic can carry 
several meanings. At a minimum, it means that there exists an unique stationary 
distribution for the process. Under this condition, the corresponding stationary 
process is ergodic in the sense of triviality of its tail cr-algebra. A second, stronger 
meaning of ergodicity for Markov processes is that for all initial distributions, the 
distribution of the process at time t converges to the (unique) stationary distribution 
as t — > 00. 
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One approach to the first kind of ergodicity involves using the stochastic equation 
to construct a "coupling form the past." Following an idea of Kendall and M0ller 
(2000), for r] 1 C rj 2 , define 

A^r? 1 ,^ 2 ) = sup X(x,r]) Ux,ri 1 ,r} 2 )= inf X(x,rj). 
Note that for ry 1 C rj 2 

\X(x, V \ V 2 )-X(x,v\v 2 )\< [ a(x,y)\ V ^ V 2 \(dy). 

We assume that N is defined on S x (— oo, oo) x [0, oo) 2 , that is, for all positive 
and negative time, and consider a system starting from time —T, that is, for t > —T 

vl ,T (B) = / Lo w^ i.T a.Tjj («)!((_,,,«,) (r)JV"(da;,ds,dr,d«) 

JBx[-r,t]x[0,=o)2 1 '- v ' ls - " 

+ / ^t+T,oo){r)V-T( dx ' dr ) 
JBx[0,oo) 

Vt T {B) = / X w i,t a .T )] (u)\t-a,co)(r)N(dx,da,dr,du) 

■/Bx[0,t]x[0,oo) 2 1 ' v ;J 

+ / l( t+ T,oo)WC-T(^,^) ; (3-1) 

JSx[0,oo) 

where we require ryi'^ C r? 2 ^. Suppose X(x,T]) < A(x) for all rj and 

c k {x)A(x)(3(dx) < oo, ft = 1,2,..., (3.2) 

which implies Condition 2.1, and suppose Condition 2.2 holds. Then we can obtain 
a solution of (3.1) by iterating 

vl" T " n+1 (B) = / \ x( i, T ,n 2,T, n)) (u)\ t - 8tX )(r)N(dx,ds,dr,du) 

JSx[-T,t]x[0,oo) 2 1 '- v s ~ s ~ 

+ / \t+T,oo){r)rilT{dx,dr) 

JBx[0,oo) 

Vt' T ' n+1 (B) = / \ J( i, T ,n 2, T , n)] {u)\ t _ ai0o) (r)N(dx,ds,dr,du) 

JBx[-T,t]x[0,oo)2 ' ' s - ;J 

+ / \t+T,oo)(r)if_:^(dx,dr), (3.3) 

,/Bx[0,oo) 

where we take r]l' T ' 1 = and 

r k ,T,1 { B ) = / \vMx)]( u )\t-s,oo){r)N{dx 1 ds,dr 1 du) 

JBx[-T,(]x[0,oo) 2 

+ / l(t+T,oo)(0^-T(*E,dr). 

JBx[0,oo) 

Note that r] 1 ' 7 '" c r) 2 ' T ' n , {77 1,T '™} is monotone increasing, and {r/ 2 ' T ' n } is monotone 
decreasing, and the limit, which must exist, will be a solution of (3.1). 

For Ccl, define (C + t) = {(s + t) : s G C}, and define the time-shift of TV 
by i? t iV(B x C x D x E 1 ) = N(B x (C + t) x D x E). Taking T = oo in (3.3), the 
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iterates 

vl'°°' n+1 (B) = / l 10iA(x ^i L oc,n )?; 2 L oo,n ) (u)l <t _ aiOo) (r)JV(da;,ds,dr,du) 

JBx(-oo,t]x[0,oo) 2 3 s 

(3.4) 

Vt'°°' n+1 ( B ) = / lf Mxv 1 ' 00 '" r,^.")](u)\t- s ,oc)(r)N(dx,ds,dr,du), 

iBx(-oo,t]x[0,oo) 2 ' ' ^ 'Is- > Is- 'I 

satisfy = H m ' n (R t N), m = 1, 2, for deterministic mappings ff m '™ : 7V(5 x 

(-oo, oo) x [0, oo) 2 ) -> and the limits 77™'°° satisfy 

riT'°° = H m (R t N), (3.5) 

where H rn = lim^oo H m n . It follows that (r)l'°°, Vt '°°) 1S stationary and ergodic. 

Lemma 3.1. Suppose that A satisfies (3.2) and Condition 2.2. Then 

l,oo v l,oo,n j 2,oo i. 2,oo,n 

Vt = lim »?t' and n t = hm r/ t ' 

ezisi and are stationary. 

Applying Theorem 2.8, any stationary solution of the martingale problem can be 
represented as a weak solution 77 of the stochastic equation on the doubly infinite 
time interval and hence coupled to versions of n l -°°- n and rj 2 -°°- n so that % ' C 
■q t C nf'°°' n , —00 < t < 00. Consequently, we have the following. 

Lemma 3.2. Suppose that A satisfies (3.2) and Condition 2.2. If 

lim f c k (x)\tf°°- n -r£-°°' n \{dx)=0 a.s. 

for k — 1,2,..., then 77 = 77 2 ' 00 = jy 1 * 00 a.s. is a stationary solution of (2.28) and 
the distribution of rj^ 00 is the unique stationary distribution for A. 

Theorem 3.3. Let A : S x Af{S) — > [0, 00) satisfy the conditions of Theorem 2.13 
with M < 1. T/ien 77 = n 2 '°° — 77 1 ' 00 a.s. is a stationary solution of (2.28) and the 
distribution of n 2 ' 00 is the unique stationary distribution for A. 

Proof. As in the proof of Theorem 2.13, 

sup C (x)i?[ / a(x,y)\ V 2 '°°' n+1 - n]^ n+1 \(dy)] 

x JS 

<M [ supc(x)E[[ a(x lV )\n 2 s ^ n -nl^idy^-^ds 

J -oo x Js 

= Msnpc(x)E[ f a(x,y)\ V 2 t '°°' n - vl'°°' n \(dy)l 

x Js 



IS 

~2,oo,n r.~A ~l,oo,n 



where the equality follows by the stationarity of 77 2 '°°> n and n 1 * 00 ^. Since the 
expression on the left is nonincreasing, its limit p exists, and we have < p < Mp. 
But M < 1, so p = 0. □ 

Definition 3.4. X(x, •) is nondecreasing, if rji C 772 implies A(x,?7i) < \(x,ri2)- 

Note that if A is nondecreasing, then for 771 C 772, A(x, 771, 772) = A(a;,77 2 ) and 
A (ai, 771, 772) = A(x,?7i). The following lemma is immediate. 



19 



Lemma 3.5. Let A be nondecreasing and satisfy (3.2) and Condition 2.2. Then 
n\'°° = liniji^oo r)l'°°' n and n 2 ' 00 = lim„_ >00 r]' 2 ' 00 ' n are, respectively, the minimal 
and maximal stationary solutions of the martingale problem for A. 

For A nondecreasing, the minimal stationary distribution can also easily be ob- 
tained as a temporal limit. 

Lemma 3.6. If uniqueness holds for (2.28) and X(x, •) is nondecreasing, then the 
process n t is attractive, that is 

77o C nl implies r}\ C rj 2 (3.6) 

for all t > 0. 

Proof. The conclusion is immediate from coupling the two processes using the 
same underlying Poisson random measure. □ 

Theorem 3.7. Suppose A satisfies (3.2) and Condition 2.2. If A(x, •) is nonde- 
creasing and 770 = 0, then the distribution ofr\t converges to the minimal stationary 
distribution. 



Proof. Note that, if we set rf_ t = 0, then rjt has the same distribution as rf , and 
by Lemma 3.6, rf s C rr\ ,co for s > —t. Since for each s > —t, rf s is monotone 
increasing in t, fj s = lim^oo 77* exists and must be a stationary process. Since 77 1 ' 00 
is the minimal stationary process, we must have fj s — nl' 00 . □ 

The same argument gives the following result on the maximal stationary distri- 
bution. 

Theorem 3.8. Suppose A satisfies (3.2) and Condition 2.2. If A(x, •) is nonde- 
creasing and 

Vo(B) = / l [0! A {x ) ] (u)1 {t _ SiOO )(r)N(dx,ds,dr,du), 

JBx(-oo,0]x[0,oo) 2 

then the distribution of r) t converges to the maximal stationary distribution. 



Remark 3.9. Note that rjo is a Poisson random measure with mean measure 
»(B) = f B A(x)(3(dx). 

We can also use the stochastic equation and estimates similar to those used in 
the proof of uniqueness to give conditions for ergodicity in the sense of convergence 
as t — > 00 for all initial distributions. 

Theorem 3.10. Let A : SxJV(S) — > [0, 00) satisfy the conditions of Theorem 2.13 
with M < 1. Then the process obtained as a solution of the system of stochastic 
equations (2.28) is temporally ergodic and the rate of convergence is exponential. 

Proof. Suppose r/ 1 and n 2 are solutions of the system (2.28) with distinct initial 
configurations t]q and t/q (equivalently, tJq and »Jq). Then, by exactly the same 
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argument as used for the proof of Lemma 2.16 we obtain 

supc(x)E[[ a(x,y)\rt-rg\(dy)] (3.7) 

x JS 

< e-*supc(a;)£;[ / a(x,y)\r& -rfi\(dy)] 

x JS 

+M [ supc(x)E[[ aix^^l-^Kdy^e-^ds. 

JO x Js 

Multiply (3.7) by e* and apply Gronwall's inequality to obtain the exponential 
rate of convergence. □ 

3.2. Spatial ergodicity. In this section, we take S = R d and assume that A is trans- 
lation invariant in the following sense. For arbitrary x,y G M. d and B G B(R d ), 
write 

T x y = x + y and T X B = B + x = {y + x;y G B}. 
Then, T x induces a transformation S x on A/"(M ) through the equation 

(S xV )(B)= V (T x B), r] e J\f(R d ),B G B(R d ). (3.8) 

Note that if rj = ^2 5 Xi , then S x rj = ^2 5 Xi - x . 

Definition 3.11. We say that A is translation invariant if A(a: + y,r/) = X(x, S y rj) 
for x,y e R d ,r] G Af(R d ). 

Definition 3.12. An A/"(M d )-valued random variable r] is translation invariant if 
the distribution of S y n does not depend on y. A probability distribution /i G 
7 , (A/'(R' i )) is translation invariant if J f(r])n(dr]) = J f (S y n) ^(drj) , for all y G R d 
and all bounded, measurable functions /. 

Definition 3.13. Let v be a translation invariant, N (R d ) -valued random variable. 
A measurable subset G C Af(R d ) is almost surely translation invariant for T], if 

1g(v) = 1g(S x ti) a.s. 

for every x G R d . rj is spatially ergodic if P{n G G} is or 1 for each almost surely 
translation invariant G C Af(M. d ). 

Similarly, for x G K d , we define S X N so that the spatial coordinate of each point 
is shifted by —x. Almost sure translation invariance of a set G G J\f(M. d x [0, oo) 3 ) 
and spatial ergodicity are defined analogously to Definition 3.13. Spatial ergodicity 
for N follows from its independence properties. 

Lemma 3.14. Suppose A is translation invariant. If r]o is translation invariant 
and spatially ergodic and the solution of (2.28) is unique, then for each t > ; r\t is 
translation invariant and spatially ergodic. 

Proof. {S y i]t,t > 0} is the solution of (2.28) with 770 replaced by Syijo and N 
replaced by S y N. By uniqueness, S y i] t must have the same distribution as i] t 
giving the stationarity. Also, by uniqueness, for measurable G C JV(M. d ) there 
exists a measurable G C Af(R d ) x J\f(R d x [0, oo) 3 ) such that 

hs^nea} = 1 {(s x ^ ,s x jv)gg} a - s - 

for all x G M. d . Consequently, spatial ergodicity for r\ t follows from the spatial 
ergodicity of (r/ ,N). □ 
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Remark 3.15. If n is temporally crgodic and 7r is the unique stationary distribu- 
tion, then it must be translation invariant since {rj t } stationary (in time) implies 
{S x rjt} is stationary. 

Lemma 3.16. Suppose that A is translation invariant and satisfies (3.2) and Con- 
dition 2.2. Then for each t, n]' 00 = lim„^oo r} t '°°' n and r] t '°° = lim„_ ) . 00 r) t '°°' n are 
spatially ergodic. 

Proof. As in (3.5), r^ 00 = H 1 (R t N) can be written as a deterministic transfor- 
mation F(t, N) of N and that Syr/l' 00 — F(t, S V N). The spatial ergodicity of r)l'°° 
then follows from the spatial ergodicity of N. □ 

Corollary 3.17. If in addition to the conditions of Lemma 3.16, A satisfies the 
conditions of Lemma 3.2, then the unique stationary distribution is spatially ergodic. 
In particular, if X satisfies the conditions of Theorem 2.13 with M < 1, then the 
unique stationary distribution is spatially ergodic. 

Corollary 3.18. If in addition to the conditions of Lemma 3.16, A is nondecreas- 
ing, then the minimal and maximal stationary distributions are spatially ergodic. 
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